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Abstract 

With a model independent method the expansion history H(z), the deceleration parameter q(z) 
of the universe and the equation of state w(z) for the dark energy are reconstructed directly from 
the 192 Sne la data points, which contain the new ESSENCE Sne la data and the high redshift 
Sne la data. We find that the evolving properties of q(z) and w(z) reconstructed from the 192 
Sne la data seem to be weaker than that obtained from the Gold set, but stronger than that from 
the SNLS set. With a combination of the 192 Sne la and BAO data, a tight constraint on ft, m o is 
obtained. At the la confidence level fl m o = 0.2781qq23> which is highly consistent with that from 
the Gold+BAO and SNLS+BAO. 

PACS numbers: 98.80.Es, 98.80.-k 
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I. INTRODUCTION 



In order to explain the present cosmic accelerating expansion discovered firstly from the 



3, H, Q, 



10] for recent reviews) 



Type la Supernovae (Sne la) [l|, |2j, , dark energy (see 

is usually assumed to exist in the universe. Dark energy is an exotic energy component with 
negative pressure, and presumably began to dominate the evolution of the universe only 
recently. Although it has been studied for nearly one decade, its nature is still puzzling. 

In general there are two kinds of methods to reconstruct the properties of dark energy from 
the observation data directly. One is to assume an arbitrary parametrization for the equation 
of state of dark energy, w(z), the potential of dark energy, V(z), the Hubble parameter, 
H(z), or the luminosity distance d^(z) with some arbitrary constants. By determining these 
constants from the observational data, we can obtain the evolving; pro per ties of dark energy. 



Different ways of parametrization have been discussed in Refs. 
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121 . 1 1 31 ] . The other is 



the no-parametric method, which usually involves directly smoothing either d^, or some 
other quantity with some characteristic smoot hing scale. Currently there are many different 
models of implementing this approach jjj], 16 ] 



Recent 
Refs. 15] 



v based on smoothing the noise of supernova data over redshift, the authors in 



16] suggested a no-parametric method in a model independent manner to recon- 



struct the expansion history of our universe and the evolving properties of dark energy. In 
Ref. [16( two kinds of supernova data: 182 Gold dataset and 115 SNLS dataset are used 
firstly to reconstruct the Hubble parameter h(z) {h(z) = H(z)/Hq) and the deceleration 
parameter q(z). It was found that both data sets give g(0) < 0, which means the universe 
is undergoing an accelerating expansion, while the Gold set seems to favor a later entering 
of this accelerating era than the SNLS one. However Gold and SNLS give a good consistent 
constraint on the present matter density parameter Q mQ (f2 m0 ~ 0.276 ± 0.023) when com- 
bined with the baryonic acoustic oscillation peak obtained from the large scale correlation 
function of luminous red galaxy in the Sloan Digital Sky Survey (SDSS) [17]. The w(z) was 
also discussed with both kinds of SNe la datasets, and it was found, in agreement with that 



obtained in Refs. 



19 



201 ] using some parameterized models, the Gold slightly favors 



an dynamically evolving dark energy with a crossing of phantom divide line while the SNLS 
does not. However, in this method, the present value of Hubble parameter, H , is needed 
prior or should be marginalized over. Since the value of Hq from different observation data 
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seems to be inconsistent and doing the marginalization wastes the compute resources, in 
this paper, we firstly generalize this model independent approach to eliminate the impact 
of Hq and then reconstruct the cosmic expansion history, H(z), the deceleration parameter, 
q(z), and the equation of state for dark energy w(z) from the new ESSENCE Sne la data. 
Beside the 162 data points given in table 9 in Ref. [6j|, which contains 60 ESSENCE Sne la, 
57 SNLS Sne la and 45 nearby Sne la, we add 30 Sne la detected at 0.216 < z < 1.755 by 
the Hubble Space Telescope |4|. 



II. THE METHOD 



Following a well known procedure in the analysis of large scale structure, Shafieloo et al. 



m 



16] use a Gaussian smoothing function rather than the top hat smoothing function to 



smooth the noise of the Sne la data directly. In order to obtain the important information of 
interested cosmological parameters expediently, In cIl(z) rather than the luminosity distance 
cIl(z) or distance module /jl(z) is studied by the following iterative method 

in 2 (m 



\nd L (zy n = lnd L (s)JU + N(z) ^(lndf (*) - lndz(*)JUi) exp 

i 

with a normalization parameter 

iV^-^^exp 



2A 2 



1) 



In 2 (i±*) 



2A 5 



(2) 



In Eqs.(l,2) A is a quantity needed to be given prior. Since a large value of A leads to a 
smooth result but depresses the accuracy of reconstruction, and inversely for a small value 
of A. So it is important to choose a reasonable value of A. Here, as in Ref. 16j, we choose 
A = 0.6. In Eq.(l), d,L{z) s n represents the smoothed luminosity distance at any redshift z 
after n iteration. When n = 1 (Il(z)q denotes a guess background model and it has been 
shown that the results are not sensitive to the chosen vaule of A and the assumed initial 

—0.9 and Q m o 



0.28 



guess model [16J. In this paper we use a wCDM model with w 
as this guess background model. \nd^ s (zi) is the observed one from the Sne la and can be 
expressed as: 

In 10 



\nd° T bs \ 



Zi) 



[H 0bs {zi) - 42.38] + ln/i = \nf obs (z t ) + \nh 



(3) 
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Here h = Hq/100 and /i o6s is the observed distance module of Sne la. Apparently using the 
above method the nuisance parameter h needs to be given prior or marginalized over. Now 
we generalize this method to eliminate the impact of h. Substituting Eq. (3) into Eq. (1), 
we obtain that 



\nd L (z) s n = In + N(z) £(ln /<**(*) ~ exp 

i 

If defining lnc^z)* = lnf(z)^ + \nh, it is easy to see 

\nf(z) s n = \nf(zy n _ 1 +N(z)J2(^f obs (z i ) - In /(z)^_ x ) exp 

i 

When n = 1 

\nf{z)\ = \nf(zr + N(z)Y,(^f obs (^-^f(zy )exp 

i 

= \nd L (z) s + N{z) J2(^f° bs (zi) ~ &d L (z) s )ex V 



2A 2 



+ In h . (4) 



In 2 



2A 5 



(5) 



^ 2 (m 

2A 2 
2A 2 



(6) 



Here dAz)^ is the luminosity distance of the suggested background model. Different from 



Ref. 



1 61 ] to iterate through Eq.(l), we use Eq. (5) to obtain the smoothed results. The 



advantage of doing so is that result is independent of h (or H ). In order to determine 
whether we obtain a best fit model after some iteration, we calculate, after each iteration, 



X 



2. 



Xr 



E 



a 



(?) 



Once this Xn reaches its minimum value we stop the iterative process and get the best fit 
result. 

By differentiating the smoothed luminosity distance we can find the Hubble parameter, 
H(z), (not h(z)) 



H(z) 



d / 100/(2) 
dz \ l + z 



(8) 



which contains the information of Hq. Then the deceleration parameter q(z) of the universe 
and the equation of state w(z) of dark energy can be obtained: 

H'(z) , 



w(z) 



q (z) = (l + z) H{z) 

[2(1 + z)/3]H'/H - 1 
\-{H Q /Hfn m AA + zf 



(9) 
(10) 
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III. THE RESULTS 



The ESSENCE program (Equation of State: Supernovae trace Cosmic Expansionan 
NOAO Survey Program) is designed to measure the history of cosmic expansion over the 
past 5 billion years. The four year data was released in Ref. 6j, which contains 60 Sne 
la points. Here we use the 162 data points given in table 9 in Ref. which contains 60 
ESSENCE Sne la, 57 SNLS Sne la and 45 nearby Sne la. In addition, as in Ref. [211 ]. we 
add 30 Sne la detected at 0.216 < z < 1.755 by the Hubble Space Telescope j^]. 

Using these 192 Sne la data points, we find when n = 42 a minimum value of \ 2 is 
obtained which can be seen from the Fig. (1). In Fig. (2) we show the reconstructed 
result of the Hubble parameter H(z) with the likelihood within la. The red line is the 
best fit result and when z = the best fit value of Ho is Ho = 65.5. Fig. (3) shows the 
evolving curves of reconstructed q(z) with la error bar. It is easy to see that the universe 
is undergoing an accelerating expansion since the present value of q(z) is less than zero, 
and the phase transition from deceleration to acceleration occurs at redshift z ~ 0.55 — 0.73 
within la, which is slightly later than that obtained from Gold set (z ~ 0.38 — 0.48) but 
earlier than that from SNLS set (z > 0.7) 16]. 

Fig. (4) shows the constraint on the present matter density parameter Q m o with Hq = 65.5 
by combining the ESSENCE Sne la and baryonic acoustic oscillation(BAO) peak obtained 
from the large scale correlation function of luminous red galaxy in the Sloan Digital Sky 
Survey (SDSS). For BAO data we use a model-independent dimensionless parameter A 
defined as 

r 1 r*i a, 1 2/3 



.4 



'mO 



/l(zi)V3 



dz 

zi Jo h(z) 



(11) 



or a flat universe, where z\ = 0.35 and A is measured to be A = 0.469(^|g) 35 ± 0.017 
17| . Here n is the spectral index of the primordial power spectrum and the WMAP3 gives 



n 



0.951 



22]. Clearly the ESSENCE Sne la and BAO give a strong constraint on Q 



mO • 



consistent with that 



At the la confidence level we obtain f2 m0 = 0.2781q 023> which is highl 
obtained from Gold+BAO and SNLS+BAO ( fi m0 = 0.276 ± 0.023) [id]. 

In Fig. (5) we plot the evolving behavior of w(z) with a marginalization of f2 m0 over 



Q m0 = 0.278+^23- The best fit (red) line shows that the ESSENCE data slightly favors 



an evolving dark energy with a crossing of phantom divide line at the near past, however 
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this evolving property is weaker than that obtained from Gold data but stronger than that 



from SNSL set obtained in Ref. 



16] . In addition, from Figs. (3) and (5) we find that the 



stringent constraint on w(z) ando(z) happens around redshift z ~ 0.5, which is consistent 



with that obtained in Refs. 



20 



23] with some parameterized models. 



IV. CONCLUSION 



In this paper, with a model independent method we have reconstructed the cosmic ex- 
pansion history and the properties of dark energy from recent ESSENCE Sne la data. We 
firstly obtain the evolution of H(z). Then the cosmic deceleration parameter q(z) and the 
equation of state w(z) of dark energy are reconstructed. Our results show that their evolu- 
tionary property reconstructed from ESSENCE data is weaker than that from Gold set, but 
is stronger than that from the SNSL one. Combining the ESSENCE Sne la and the BAO 
data, a tight constraint on f2 m o is obtained. At a la confidence level fi m o = 0.2781q;q23, 
which is highly consistent with that obtained from Gold+BAO and SNLS+BAO. Remark- 
ably, as that obtained with some parameterized model {20 . 
and q(z) seem to happen at about z ~ 0.5. 



231 ] . the tight constraints on w(z) 
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FIG. 1: Computed x 2 f° r the reconstructed results at each iteration for the 192 Sne la data. 
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FIG. 2: The reconstructed evolutionary curves of the Hubble parameter H{z) with the likelihood 
within la. The red line is the best recovered result. 
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FIG. 3: The reconstructed evolutionary curves of the deceleration parameter q(z) with the likeli- 
hood within la. The red line is the best recovered result. 
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FIG. 4: The constraint on f2 m o from the combination of 192 Sne la and BAO data. The red and 
black lines show the derived value of A/\/Q m Q from the 192 Sne la dataset within la and the blue 
lines are the results from the BAO data. 
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